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We use information from rare nonleptonic decays of heavy-quark mesons to put new 
bounds on the magnitudes of certain product combinations of baryon nonconserving R- 
parity violating couplings in supersymmetric models. Product combinations of lepton 
and baryon nonconserving R-parity violating couplings are also considered in the light of 
existing bounds on nucleon decay. Contrary to popular impression, a few such combina- 
tions are shown to remain essentially unconstrained. 
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Though the minimal supersymmetric standard model (MSSM) JTJ] is a leading candi- 
date for new physics beyond the standard model, the conservation of R-parity, R p , which 
is assumed in the model has no real theoretical justification. This has motivated many 
authors to consider alternatives in which R p is explicitly broken. In such models, 
sparticles can decay into non-supersymmetric particles alone, leading to novel signatures 
in search experiments and unusual decay processes. 

The most general i? p -violating superpotential that one can write with the MSSM 
superfields, in the usual notation, is 

W = A, ; /,/., /.,/-.'/, + X'ijkLiQjDk + X'^j;/),!), (1) 

Here, i, j, k are generation indices and we have rotated away a term of the form fi^LiHy 
Since the A^ term is symmetric under exchange of i and j, and antisymmetric in color, 
it must be antisymmetric in flavor, thus we have A^ = —Xjik- Similarly, X'-j k = —\ikj- 
The number of couplings is then 36 lepton nonconserving couplings (9 of the A type and 
27 of the A' type) and 9 baryon nonconserving couplings (all of the A" type) in total. 

It is generally thought that A, A' type couplings cannot coexist with A" type couplings 
since both baryon and lepton number violations would lead to too rapid a proton decay. 
For this reason, previous authors have considered either lepton nonconserving or baryon 
nonconserving couplings, but not both. We will first make this assumption and consider 
the baryon number violating A" couplings alone, since most of the earlier effort has 
been focused on A and A' couplings ||. Later, we will examine proton decay in the 
presence of all three types of couplings. Severe constraints on violating couplings can 
be obtained by requiring that the cosmological baryon asymmetry not be washed out|Q, 
but it is possible to evade these bounds p. 

Our philosophy throughout will be that we expect the couplings involving third gen- 
eration fields to be the largest. There are two reasons for this. The only other Yukawa 
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couplings in the model (the Yukawa couplings of the quarks to the Higgs fields) exhibit a 
very strongly hierarchical generational structure and thus one would expect the R-parity 
violating couplings to do so as well. Second, many models result in the scalar top being 
the lightest of the squarks, and thus processes involving the scalar top may not be as 
suppressed by large propagators. We will thus concentrate on couplings involving higher 
generations, but will keep our bounds as general as possible. 

We can write the 9 different A" couplings as A£ 6s , X'{ bd , \'{ sd , X'^ s , X% d , X'^ sd , X'^ bs , X'l bd and 
X'usd- Let us first recount the existing constraints on these couplings. Brahmachari and 
Roy |J showed that the requirement of perturbative unification typically places a bound 
of between 1.10 and 1.25 on many of the couplings. This was generalized to all the cou- 
plings by Goity and Sher ||. The latter also showed, following earlier work || ||, that 
\X" sd \ can be strongly bounded by the nonobservation of double nucleon decay into two 
kaons (such as 16 — > 14 C K + K + , which would have been seen in water Cerenkov detec- 
tors), and |A" M | can be strongly bounded by the nonobservation of neutron-antineutron 
oscillations. Their bounds, for squark masses of 300 GeV, were \X" bd \ < 5 x 10 -3 and 
\X'u Sd \ < 10~ 6 . In the work of Barbieri and Masiero |J, some bounds on products of 
couplings were obtained by considering K-K mixing; these bounds will be discussed 
shortly. Finally, bounds from the bb induced vertex correction to the decay of the Z 
into two charged leptons have recently been obtained [10]; though potentially interest- 



ing, with present data they are not significantly better than the bound from perturbative 
unification. 

In this Letter, we note that many additional and interesting bounds on the A" cou- 
plings can be obtained by considering rare two-body nonleptonic decays of heavy-quark 
mesons. We shall mostly consider B decays, but also, in some cases, D decays. Let 
us begin by considering the implications for A" couplings from such processes.. Since 
lepton number is assumed to be conserved, only AB = and AB = 2 decays can occur. 



Thus any bounds will be on the products of two A" couplings. Furthermore, since any 
B-decay will change the number of "b-flavors" by one unit, bounds from there will ap- 
ply to products of the form X Ui b s X Uj sd or X Ui bdX Uj sd- We first consider AB = (baryon 
number conserving) decays which actually give the best bounds and later comment on 
the AB = 2 processes. 

In our calculation of two-body nonleptonic decays of heavy-quark mesons, we follow 
the computational method of Carlson and Milana |Tl[ which is based upon the formalism 



of Brodsky and Lepage | 12| . First, we neglect all light meson masses (7r's, K's). Then we 



make use of the fact that the relative momentum between the quark and the antiquark 
for each of the qq pairs within the decaying meson and within each of the final state 
mesons is low. The large quark momentum transfers, needed to move a quark from one 
meson to another or to produce a qq pair which enter different mesons, can be caused 
either by single-gluon exchange or through the emission of a virtual squark. However, 
both are needed in each diagram for the correct distribution of relative momentum, so 
we consider diagrams that involve both. Alternatively, one could consider the full meson 
wave functions ]13[ with the tail of the wave functions being crucial; if properly done, this 
should be equivalent to the previous calculation since at high enough relative momentum 
the tail is indeed given by 1-gluon exchange. 

First, take the decay B + — > K°K + (or equivalently, B~ — > K°K~). This has an 
extremely small rate in the Standard Model, being penguin- suppressed and also reduced 
by the small CKM element V u b in the amplitude. The dominant diagrams contributing 
to this process for nonzero A" couplings are shown in Fig. 1. In each of these diagrams, 
the gluon is spacelike and the squark (of charge 2/3) is timelike. This generates from the 
gluon propagator an overall enhancement factor of mg/ (m# — 7715) ~ 10 in the amplitude, 
mg and mj, being the B-meson and b-quark masses respectively. (This factor is just the 
inverse of the fraction of the B- momentum assigned to the light quark.) One can draw 



similar diagrams interchanging the squark and gluon internal lines and appropriately 
relabeling the quark lines. Each of these latter diagrams would have a timelike gluon 
and a spacelike squark (now of charge —1/3). For these, however, the overall factor 
of mB/{iriB — m&) does not materialize so that contributions from these diagrams are 
significantly subdominant and can be neglected. It is of interest to point out that the 
decay B° — > K + K~~ can proceed only through these latter diagrams, with the squark 
having charge 2/3, and therefore yields a significantly weaker bound on the same product 
of A" couplings than does B + — > K°K + , in spite of stronger experimental limits on the 
branching ratio. 

Our result is most conveniently expressed as a ratio of T(B + — > K°K + ) to the partial 
width of another B + decay channel (specifically B + — > K + J/ip) that proceeds unsup- 
pressed in the Standard Model. This description eliminates many of the uncertainties in 
the coefficient factors. We find, considering only t-squark contributions, that 

K°K+) _ ( m% A - 1 / f K y \K b Xs d \\rn w / mi y 
T(B+^K+J/ii) \ L ml) [fj/J (G F m* w y\V cb ms\ 2 [ } K) 

Here, mi is the mass of the scalar top, mju is the mass of the J/ip meson, fx and /ju 
are the decay constants of the K and J/ip, which are related to their wave functions at 
the origin, and V cb and V cs are CKM elements. We shall use jji^j fx — 2.55. Using the 
experimental branching ratio for B + — > K + J/tp, namely 10.2 x 10~ 4 , we have 



B.R.(B + - K°K+) ~ 1.97\\'; bs \'; sd \ 2 (m w / mi r. (3) 
On using the recent experimental upper bound [15] of 5 x 10~ 5 on the branching ratio, 



and noting that one can replace the scalar top with a scalar charm or scalar up, we then 
have 

IV V < 5 x 10- 3 , (4) 
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for q = t,c, u. We have redone the computation of Eq. 2 using the methods of heavy 
quark symmetry and find a 15% downward revision in the upper boundQ. 

One can repeat the calculation for the decay B + — ► K°tt + (or equivalently B~ — ► 
K°7r~) in much the same way. The result (with only the t-squark contribution being 
considered) is 

B.R.(B+ - K%+) = l.^\\l d \'U 2 {rn w /m~) A . (5) 

Again, using the experimental upper bound [15] of 5 x 10~ 5 on the branching ratio, we 
have 

\\ bd \ sd \m w _ < 4 A x 10 _ 3; (6) 



for q = t,c, u. 

Though these methods of calculation would be less reliable for two-body nonleptonic 
decays of D-mesons, we can use a similar approach there to get an order of magnitude 
estimate of the corresponding bounds. However, we find that the bounds are significantly 
higher than bounds obtained from D — D mixing, discussed below. 

Additional bounds were obtained by Barbieri and Masiero || from the contribution 
of K-K mixing to the K^-Ks mass difference. There are two main box diagrams towards 
this contribution, shown in Fig. 2. Assuming that one or the other of these diagrams is 
dominant, the results of Barbieri and Masiero lead to | 

\K b Xu\ < mm (6 x 10-\m- t /m w ), 3 x 10^(m- t /m w ) 2 ) (7) 

\KbsKJ < mm (6 x 10- 4 K/m w ), 2 x lQ-\m~ c /m w f) (8) 

These authors did assume that the top quark was much lighter than scalar top, and that 
^^This gives some idea of the theoretical uncertainty in the bound. 

2 We are considering bounds on the real part of the coupling constants only. Bounds on possible 
imaginary parts were also discussed by Barbieri and Masiero. If the couplings have significant imaginary 
parts, then in the bounds in this paper arising from K-K and D-D mixing, |A a Af,| must be replaced by 
|Re(A: 2 A2|)|V2. 
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all squark masses are degenerate. The former assumption has been invalidated by recent 
data, so their bounds need to be revised. 

These results can be generalized by noting that the exchanged quarks in Fig. 2a can 
be any two charge-2/3 quarks, and also that the exchanged quarks in Fig. 2b can consist 
of one c-quark and one t-quark (w-quark contributions are suppressed by a mass insertion 
in this diagram). We have calculated the possible contributions, assuming that all of the 
squark masses are equal except that of the scalar top, and using a mass of 175 GeV 
for the top quark with updated CKM angles.^ In Fig 2a, we have also included (as did 
Barbieri and Masiero ) the contribution arising by replacing all of the particles in the box 
with their superpartners; we have not included the similar contribution from Fig. 2b due 
to the extra unknown parameters (it is unlikely that this contribution will almost exactly 
cancel the calculated contributions; thus our results will not be significantly affected by 
them). We require that these contributions not exceed the standard model contribution 
(which is uncertain by roughly a factor of two), and have plotted the upper bound for 
various products of couplings in Fig. 3. It is not hard to see that the contribution from 
B-B mixing to the B^-Bs mass difference will give bounds on the same couplings, but 
will be weaker. 

Bounds from D-D mixing can also be considered. We find that one of the two box 
diagrams is suppressed by small CKM angles and the other gives 

\K bs K bs \<3.1xlO- 3 (m 5 /m w ) (9) 

What about AB = 2 decays? One can envision the process of Figure 4, which will lead 
to the decay B — > or A A . However, a simple estimate of the rate gives branching 

ratios (assuming the B-violating couplings are unity and the scalar quark masses are 

near the W mass) of O(10~ 8 ). The smallness of the rate is due in large part to two 

3 In the limit of small top quark mass, our analytic result for the effective Hamiltonian is a factor of 
two smaller than that given by Barbieri and Masiero. 
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small CKM elements in the amplitude. Thus, such processes will not provide interesting 
bounds unless 10 10 B-decays are studied. 

It is generally assumed that the presence of both lepton nonconserving terms and 
baryon nonconserving terms leads to unacceptably rapid proton decay. However, if 
enough third generation fields are involved, proton decay can be sufficiently suppressed 
as to make some of the bounds very weak (or, in a few cases, nonexistent). To see this, 
suppose both A and A" terms exist. Consider the bound on the product |A MTT A^ sd |. This 
will lead to proton decay through the diagram of Figure 5a. Although there is a sup- 
pression due to mixing angles and heavy squark propagators, the proton lifetime bound 
gives a strong bound of |A MTr A" sd | < 1CT 14 . This bound is independent of the final state 
leptons, and thus applies to all 9 of the A couplings. Similar bounds can be obtained for 
all A" couplings with at most one heavy field (which is then the scalar quark); we obtain 
\\ ijk K bd \ < 10- 13 , \\ ijk K bs \ < 10- 12 , \X ijk KJ < 10- 13 and \\ ijk \'l sd \ < 1(T 12 . However, 
if the A" coupling has two heavy fields, a loop is necessary, as shown in Figure 5b. This 
gives much weaker bounds; we obtain \\ijkKbs\ < 10 -2 , Wjk^tbdl < 10~ 3 , \KjkKbs\ < 10 -3 
and \XijkKbd\ < 10~ 2 . We thus see that the lack of obvservation of proton decay does 
NOT always give very strong bounds on the product of the lepton number violating and 
baryon number violating couplings. 

Finally, we consider the product of A' and A" couplings. Here, there are 27 x 9 possible 
products, of the form |A^ fe A" bc |. The diagrams that can lead to proton decay are shown 
in Fig. 6. For each of these, one can have a r or c, b, t quark on an external leg, in 
which case that leg must be virtual and decay through a W. We have examined all 
posssible products of couplings and found that the vast majority are tightly bounded 
(product is less than 10~ 6 ), but some are not. Rather than list the bounds for all 243 
combinations, only the bounds which are greater than 10~ 6 (for the product of the A' and 
A" couplings) will be given explicitly. It is found that all products with a A" sd , \" bd and 
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\'l bs are smaller than 10~ 9 . The same is true for \" sd , \" bd and \" bs , exception \X' tbl X" sd \ 
which is < 10 _1 , |A^A^, d | which is < 1CT 2 and \X' udl \% s \, for which no bound better than 
the unitarity bound could be found. Here I is any lepton. For X" sd , the only bound which 
is not very small is the combination |A^A" sd |, which is bounded only by unitarity For 

Kbd an d Kbsi we fi n d \Kd(e,ti)Kbd\ < 10 2 ' \Ks( e ,iJ,)Kbs\ < io~ 2 ' \K(d,s)(e,ti)Kbd\ < io~ 3 > 

\K( d ,s )M X " b s\ < 10- 3 , \\' cdT K bd \ < 10- 5 , |A^A&.| < 10- 5 , whereas | A' ( „ )C)dr A^l and 
|AL c ) sr A^| are bounded only by unitarity. All other bounds are quite tiny. Thus, of 
the 243 combinations of couplings, thirty have bounds greater than 10 -6 , and eight are 
completely unconstrained by the lack of observation so far of proton decay. 

Could the unconstrained couplings be bounded by rare B decays? One can envision 
a diagram similar to Fig. 1 in which the upper vertex is a A' vertex; this would lead to 
B decay into three quarks and a lepton, such as B — > pr. However, in all of the above 
30 couplings, the scalar quark leaving the A' vertex is different from that entering the A" 
vertex, and thus a loop will be necessary, suppressing the rate. As originally noted by 



Thorndike and Poling[16], therefore, the bounds from proton decay will be better than 
those from B decay in all cases, unless branching ratios of 10~ 8 or better are obtained. 

It is interesting that the standard supersymmetric model can contain some baryon 
number and lepton number violating coupling constants which are of order unity, and 
which do not lead to excessively fast proton decay. Such couplings could be measured 
when squarks and sleptons are discovered, since they will lead to baryon and lepton 
number violating squark and slepton decays. These have been extensively discussed by 
Dreiner and Ross and by Dimopoulos, et al. JlTj , who analyze the impact of such decays 
on phenomenology. 

P.R. would like to acknowledge the hospitality of the College of William and Mary and 
of the CEBAF laboratory during the course of this work. We would like to thank Michael 
Peskin and Ed Thorndike for useful discussions, and Frank Wurthwein and Karen Lingel 
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Figure 1: The dominant diagrams contributing to B + — > K°K + . The t refers to the 
scalar top quark, the g to the gluon. 



Figure 2: Contribution of R-parity violating couplings to the K L -K$ mass difference. 



Figure 3: Bounds on products of A" couplings, plotted vs. the scalar bottom mass, m s , 
arising from constraints from K-K mixing. The ratio of the scalar top mass to the scalar 
bottom mass is taken to be 0.5; all other scalar quark masses are degenerate with the 
scalar bottom mass. The bounds on |A" M A" bs |, |A" bs A" bd |, lA'^A'^J are given by the solid, 
dashed, and dotted lines respectively. The bound on (A^A^A^A"^ 1 / 2 is y/2 lower 
than the dotted line; the bound on |A^A^ M A^ A LI and \KbdKbsKbd X tbs\ are S iven b Y 
the products of the solid and dotted lines, the bound on |A" bd A" bs | is also the dotted line; 
and the bound on |A" bd A" bs | is larger than that of |A" bs A" bd | by a factor of roughly 1.7. If 
the scalar top mass is taken to be 0.9m s , then only the |A" bd A" b J curve changes, increasing 
by roughly 50% over the entire range. 



Figure 4: A typical contribution of R-parity violating couplings to the baryon number 
violating decay of the B°, in this case yielding B° — > AA. 



Figure 5: Diagram (a) shows a typical proton decay arising from the existence of both A 
and A" R-parity violating terms. In diagram (b), the contribution is given for the case 
in which the A" term contains two heavy quarks. 



Figure 6: Contributions to proton decay arising from the existence of both A' and A" 
R-parity violating terms. The tree-level diagram contributes if the A' and A" terms each 
have a single, identical heavy field; in other cases, either the loop diagrams contribute or 
one of the external lines is a heavy quark which must then emit a virtual W (or, in some 
cases, both must occur). 



12 



